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Tomographic Magnification of Lyman Break Galaxies in 
The Deep Lens Survey 
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ABSTRACT 

Using about 450,000 galaxies in the Deep Lens Survey, we present a detection of the 
gravitational magnification of z > 4 Lyman Break Galaxies by massive foreground 
galaxies with 0.4 < z < 1.0, grouped by redshift. The magnification signal is detected 
at S/N greater than 20, and rigorous checks confirm that it is not contaminated by any 
galaxy sample overlap in redshift. The inferred galaxy mass profiles are consistent with 
earlier lensing analyses at lower redshift. We then explore the tomographic lens mag- 
nification signal by splitting our foreground galaxy sample into 7 redshift bins. Com- 
bining galaxy-magnification cross-correlations and galaxy angular auto-correlations, 
we develop a bias-independent estimator of the tomographic signal. As a diagnostic 
of magnification tomography, the measurement of this estimator rejects a fiat dark 
matter dominated universe at > 7.5a with a fixed as and is found to be consistent 
with the expected rcdshift-dependence of the WMAP7 ACDM cosmology. 

Key words: gravitational lensing: weak — large-scale structure of the Universe — 
cosmology: observations — galaxies: high-redshift — galaxies: haloes 



1 INTRODUCTION 

Lensing tomography is a powerful probe of cosmology. By 
comparing the amplitude of lensing effects across multiple 
redshift baselines, we can observe both the geometric expan- 
sion of the universe and the growth of structure as a function 
of time. These are two key diagnostics of dark energy, mak- 
ing tomography an essential p art of the toolkit iden tified by 
the Dark Energy Task Force |Albrecht et al. I [20061 ) for use 
in upcoming sky surveys. 

Gravitational lensing has two observational signatures: 
shear and magnification. While gravitational shear has been 
the main focus of observational efforts in the past decade, 
magnification measurements are now possible with the avail- 
ability of deeper surveys together with stable photome- 
try. Weak magnificat ion by galaxies robustly detected by 
IScranton et all (|2005l ) based on the cross-correlation be- 
tween low redshift galaxies with distant quasars from the 
Sloan Digital Sky Survey (SDSS). Subsequent work by 
I Menard et al.l (|2010l ) showed that shear and magnification 
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measurements give consistent mass profiles around galax- 
ies. Since then, the detection of gravitational magnification 
has been reported in various contexts: usin g Lyman Break 
Galaxies (LBGsl dffildebrandt et all 120091 ), sub- millimeter 



galax ies l|Wang et al.l 12011 1 ). galaxy groups (|Ford et al.l 
|2012| ), and clusters (|Hildebrandt et all 1 1201 ll ). Several es- 
tima tors have been used , based on source number den- 
sity (IScranton et al.ll2005|). brightness chang e (I Menard et al 
2010 



n et al. zuuo), brightness change (Menard et al 
change (iHuff fc Graves! l201ll: iSchmidt et al 



2012) and quasar variability ( Bauer et al.ll201lf ). 



So far attempts to detect a tomographic lensing sig- 
nal have only been m ade through shear measurements 
(jSchrabback et all l2010h . Such analyses require accurate 
shape measurements as a function of redshift. In contrast, 
magnification measurements do not rely on shape measure- 
ments and can even be done for unresolved sources. In this 
paper, we use LBGs at redshift z ~ 4 as background sources 
and cross-correlate with foreground lenses. We then extend 
the magnification measurement, breaking the foreground 
lenses into bins of photometric redshift. By selecting mul- 
tiple photometric redshift samples, we probe the growth of 
structure and distance for several different ranges of cos- 
mic time, albeit at the expense of having galaxy samples 
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with redshift-dependent biases. To account for this bias, 
we normalize each of the lensing measurements by the ob- 
served foreground sample autocorrelati on, resulting in a sig- 
nal that is nearly bias independent |jain fc Tavlorl 120031 ; 
van Waerbekeilioiol ). 

The paper is structured as follows; in Sj2]we derive the 
magnification formalism and the expected tomographic sig- 
nal; Sj3] covers the details of the Deep Lens Survey (DLS), 
the data selection criteria, and the analysis techniques used; 
SJH describes the results; and in we summarize and dis- 
cuss future directions. We discuss detailed survey parame- 
ters, possible systematics, and theoretical modeling of the 
expected signal in the Appendix. 

All cosmological calculations assume, unless otherwise 
stated, a fiat universe with SIa = 0.73, Q m — 0.27, and 
erg = 0.8. Magnitudes and color cuts are done using AB 
magnitudes. 



2 THEORY 

Gravitational magnification can change th e apparent den- 
sity of background sources (jNarava n 1989) and as a result 
induce apparent cross-correlations between foreground and 
background galaxy populations. In this section we present 
the formalism des cribing magnification-induced sp atial cor- 
relations following iBartelmann fc Schneider! (|200ll ). 

For a flux-limited sample, the integrated number of 
galaxies above a given flux threshold / can be approximated 
as 



No(> f) ~ Af~ 



(1) 



where Ao is the number of galaxies in area A and a is 
the power-law slope of the number counts. Magnification 
increases the depth observed by reducing the effective flux 
limit to f —¥ f/n and simultaneously decreasing the area 
of the survey by A — > A/fi, where y, is the magnification. 
Combining these effects, the relation for the galaxy number 
counts becomes 



N(> f)~~A[$- 



^ _1 Ao(> /) 



(2) 



From this functional form, we see that for values of a > 1 the 
number of sources observed is increased while the opposite 
is true for a < 1. At the critical value of a = 1, no effect 
from magnification is observed. 

We introduce the cross-correlation between lenses and 
sources as follows: we can define the over-density of sources 
N s as a function of position on the sky <j> to be 



Ns(<j>)-{N a ) 
(Ns) ■ 



(3) 



where (N s ) is the average density. Under the effects of mag- 
nification given in Equation [21 the number of sources is 



N s (<j>) = n(c/>r- 1 {N a }. 



(4) 



We consider small departures of fj, from unity and substi- 
tute fj, = 1 + 8/i where |<5 M | <C 1. Taylor expanding, we can 



write /i 



1 + (a — 1)<S M , thus the over-density under 



Hence, the angular cross-correlation between lenses and 
sources induced by magnification can be written as 



w ls (6) = (a-l){5»{<f>)5 g {<f> + e)) 
= (a-lK i( («) 



(6) 



where 6 is the angle from the lensing galaxy center and the 
foreground galaxy density is 8 g . We use Limber's approxima- 
tion to convert the intrinsically 3D distribution of galaxies 
to projected angle. We do this by e xpanding the above equa- 
tion in the same framework as in IBartelmann fc Schneider] 
We then have 



Wf.Sg (0) 



kdk 
~2^~ 



b(k, X )r(k, z)Pum(k, X )Mxke) (7) 



where x ls the co-moving distance, Pdm is the dark matter 
power spectrum, b is the galaxy bias as a function of scale 
and distance, r is correlation coefficient also as a function so 
scale and distance, r/i is the co-moving distance distribution 
of the foreground lensing galaxies, Jo is the zeroth order 
Bessel function, and 



K{x) 



3H fl m x 



dxVs(x 



X 



X 



(8) 



magnification becomes 



is the lensing kernel-weighted distribution of background 
sources rj s (x)- 

For the remainder of the paper, we make the assumption 
that the galaxy bias is linear, that is b(k, x) — ► b, a constant, 
and r(k, x) — ► 1. There are limitations to these assumptions. 
If the bias is stochastic, r < 1, and the bias is a strong func- 
tion of k then interpreting these measurements cosmologi- 
cally will be diffi cult. There is observational evidence tha t 
this is the case l|Hoekstra et al.l |2002| ; ISimon et ail 120071) . 
howe ver a recent measurement in COSMOS ( Julio et al.l 
l2012h found that for galaxies selected in photometric red- 
shift, at redshifts similar to this analysis, b was linear to 
small scales and r = 1 to within the measurement error. 
We attempt to mitigate these effects in the next section and 
show how well we can for realistic HODs, similar to our own. 



2.1 Tomographic Signal 

The galaxy bias b complicates the interpretation of mag- 
nification results with regard to cosmology and large scale 
structure, since the bias is dependent on galaxy type, bright- 
ness, and mass - any or all of which may vary in the lensing 
sample as a function of redshift. If we want to isolate the 
cosmological evolution of the lensing signal, we have two 
choices. Either we can model the galaxy bias for each red- 
shift bin and marginalize over the parameters for that model 
or we can try an entirely empirical approach where we at- 
tempt to cancel out the bias using the autocorrelation of 
the foreground sample by constructing a bias independent 
(or nearly bias independent) estimator. For the purposes of 
this paper, we choose the latter. 

Using the same assumption that the density of galaxies 
follows the dark matter over-density as bS, the autocorrela- 
tion for the lensing galaxies is given by 



5,(0) = (a -!)$„ 



(5) 



w u {6) =b 2 {8{4>)5(<t> + 9)). 



(9) 
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Using Limber's approximation, this becomes 



w n (6) = b dxmix) 
Jo 



krlk 



(10) 



where the definitions follow from the previous section. In 
this linear approximation of the galaxy bias, we can remove 
the bias by taking an appropriate ratio of the magnification 
and autocorrelation signals: 



Wtl 



(11) 



This ratio is independent of galaxy bias, as long as the linear 
approximation of the bias holds. 

Observationally, we do not directly measure wu , instead 
we measure the true autocorrelation with an additional con- 
stant value. Pair conservation introduces a measurement 
bias C: 



c{6) = [i + Wll (e)]N- 2 Y, 



(12) 



for N galaxies in the survey area with the sum taken over 
all pairs. In practice, C(9) is nearly constant as function 
of scale, resulting in an ov erall suppression of wu (| Peebles! 
Il98d ; IScranton et al.ll2002T ). We use Equation [T2l to correct 
the amplitudes of the measured auto-correlation in the cor- 
relation ratio, giving us a final form of 



K 



Wu 



C 



(13) 



We model Equation [13] in the quasi-linear and non- 
linear regime, and also account for the effects of the integral 
constraint on the autocorrelation measurements. For this 
purpose, we use the code NICAEjB 

to generate the dark 
matter power spectrum, Pdm, and subsequent correlations. 

To test the bias independence of the quan tity TZ at sma ll 
scales we utilize the halo model as defined in ISeliaJj l|200d) . 
To estimate the galaxy power spectrum, P gg , and galaxy- 
matter cross power spectrum, P gm we use two models for 
populating galaxies in h alos. The first is the halo occupation 
distribution (H OP) from | Zheng et al.l ((2007) with model pa- 
rameters from I Wake et al.l (1201 if). For t he second we use 
the HOD from iMandelbaum et all |2005t ) with parameters 
as measured in the DLS from galaxy-galaxy lensing (private 
communication A. Choi). We chose these parameters as they 
bracket this analysis both in redshift and halo masses. Fig- 
ure [T] shows the resultant ratio TZ at two physical distances, 
1 Mpc/fo and 0.5 Mpc/ft, for the two HODs compared to 
the dark matter. It should be noted that the redshift bin- 
nings used in this plot are not the same as those from the 
data, however, the resultant bias independence of the ratios 
will be similar. The ratio using HODs agree with the dark 
matter only ratio, at worst deviating by 20%. The TZ and TZ 
ratios are then bias-independent at the scales used in this 
analysis. For full details on the HODs and halo modeling see 
Appendix [Al 



1 Kilbinger, Martin: http://www2.iap.fr/users/kilbinge/nicaea/ 
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Figure 1. Predicted curves for the ratio TZ(z) on 0.5 and 1.0 
Mpc//i scales, for a range of models. The solid black line is the 
dark matter-only ratio, dashed curves use HODs derived from 
DLS galaxy-galaxy lensing at z ~ 0.6 (private communication 
A. Choi), an d dot-dashed curv es are fits to galaxy clustering at 
z ~ 1.1 from lWake et al.1 |201lh . For both HOD classes, the red, 
green and blue curves run over a range of galaxy bias. See Table 
lAll for more details. 



3 DATA 

Reliably measuring weak lensing magnification demands ac- 
curate, consistent photometry over the whole survey area 
utilized. In this section we lay out the calibration steps we 
use to achieve this, and discuss the selection of both lenses 
and sources. 



3.1 The Deep Lens Survey 

The Deep Lens Survey (DLS) (|Wittman et al.ll2002r ) is a four 
band, BVRz, CCD survey of 20 square degrees observed on 
both the KPNO Mayall and CTIO Blanco telescopes, opti- 
mized for the observation of gravitational lensing. The sur- 
vey is composed of five 4 deg 2 fields (two northern and three 
southern) , widely spaced in RA/Dec at high galactic lati- 
tude . The MOSAIC imagers (|Pogge et al.ll99Sl ; fMuller et al.1 
1998) served as the observing instruments during the course 
of the survey, with each field consisting of a 3 x 3 grid of 
MOSAIC footprints for a total of 45 subfields. Nights with 
the best seeing (PSF FWHM < 0.9") were reserved for R 
band, leading to a total exposure time of 18,000 seconds. The 
remaining bands all have a total exposure time of at least 
12,000 seconds. Because V observations were often taken on 
good seeing nights after finishing R observations, the seeing 
in V band is also ~ 1". The survey is 50 % complete in 
terms of object recovery to 26ab in R band, 25.5as in BV 
and 24.5ab in z. 

The DLS image processing and photometry are de- 
scribed in detail in an upcoming data release paper 
(Wittman et al. in prep.) as will be the photometric red- 
shifts used in this analysis (Thorman & Schmidt in prep.). 
The details of these pipelines relevant to this analysis can 
be found in the Appendix. 
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PRIMUS Spectroscopic Redshift z 

Figure 2. Plots of the DLS photometric redshift performance. The figure on the left shows the summed redshift posterior probability 
p(z) vs spectroscopic redshift from ~ 9000 PRIMUS survey galaxies binned by spectroscopic redshift. The contour values are linearly 
spaced and range from 200/dz 2 to 3800/dz 2 where dz is 0.01 in redshift . Integrating over the full range of spectroscopic and photometric 
redshifts gives the total number of galaxies in PRIMUS that overlap the DLS. The right panel shows the redshift distribution for the 
7 photometric redshift-selected bins based on the aggregate p(z) for the galaxies in each bin, weighted by the likelihood that a given 
galaxy is inside the bin. The curves are normalized to a maximum of one. 



3.2 Foreground Lens Selection 

To measure magnification and tomography the lens sample 
must be free of spatially varying survey systematics (see Ap- 
pendix [JC]for more details) as well as have good photometric 
redshifts. We accomplish this by selecting galaxies with mag- 
nitudes in the range 20 ^ run ^ 24, which is well below and 
well above the saturation and detection limits respectively. 
This cut also assures quality photometric redshifts due to 
the high signal to noise of the galaxies in this range. 

We used the redshift posterior probability, p(z), cal- 
culated by the Bayesian redshift estimator code, BPZ for 
each galaxy (see i]B2l for details) to define redshift bins. 
The posterior probability distribution for each galaxy is 
summed, yielding an estimate of the redshift distribution 
for the foreground lenses (albeit one that was convolved with 
the redshift-dependent scatter of the photometric redshifts) . 
We then split this distribution into 15 bins of equal likeli- 
hood over the range < z < 5. Since the summed likelihood 
is a convolution of the true redshift distribution and the 
photometric redshift scatter, this binning produces broader 
redshift bins at higher redshift, as one would expect. The 
choice of 15 bins guarantees that no bin is narrower than 
8z ~ 0.06, the expected scatter in the DLS photometric red- 
shifts. 

From the initial binning, we select 7 bins that span the 
peak of the lensing kernel for a redshift z = 4 (the redshift 
of the Lyman Break Galaxies, see M3.3p source. The filter set 
of the DLS, BVRz' , allows for accurate redshifts within the 
range z = 0.4— 1.0, as such we only consider galaxies within 
this redshift range in the analysis. 

To test how well the photometric redshifts ar e perform- 
ing w e use the spectroscopic survey PRIMUS l|Coil et al.l 
i|201ll ). Cool et al. in prep.) as a cross check. The left panel 
of Figure [2] shows the summed p(z) versus spectroscopic red- 
shifts from the PRIMUS survey that overlap the DLS. From 



PRIMUS we have ~ 9000 spectroscopic redshifts which are 
100% complete to an tur band magnitude of 22.8ab and 
30% complete to ma of 23.3yis. Contours of p(z) in the 
redshift range z = 0.4 — 1.0 track the one-to-one line in a 
mostly unbiased manner, with their expectation value within 
z = 0.02 of the mean spectroscopic redshift value. There are 
also no significant spurious peaks in range of redshifts plot- 
ted, however, there are degeneracies between redshifts of 
z ~ 2 and z < 0.3 caused by a lack of U band data. We find 
that, for the range of redshifts used in this analysis, the full 
p(z) is robust against catastrophic outliers and we find sim i- 
lar benefits to using the full p(z) to those in lWittmanl ([2009) . 
For the galaxies with spectroscopic redshifts from PRIMUS, 
the catastrophic outlier rate for point estimate redshifts de- 
fined as the peak of the redshift posterior is 7.6%, and drops 
to 3.5% in the z — 0.4 — 1.0 range used in this analysis. 

For each galaxy, we calculated the probability Pi that it 
is within a given redshift bin i's bounds z m in,i < z < z max ,i 

Pi = p(z) dz. (14) 

To select that a galaxy be in a given bin, we require that Pi > 

0. 16, which is the single tailed la probability of being in bin 

1. To avoid double-counting, we then weight each galaxy by 
their Pi value for a given bin. This selection is similar to 
a photometric redshift ODDS cut, where a requirement is 
made that the integrated probability within a given range 
of the peak of a galaxy's redshift posterior is above some 
base threshold. Defining a catastrophic outlier rate as the 
excess probability outside of the range zi ow — 0.15 * (1 -\-z) - 
Zhigh + 0.15 * (1 + z), where Zi ow and Zhigh are the lower and 
upper bin bounds respectively and z is the average redshift 
of the bin, we find that probability outside of this range is 
~ 3% for each bin which similar to that of the outlier rate 
estimated for the peak redshift. 

The right panel in Figure [2] shows the summed p(z) 
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Table 1. Redshift bin properties. 
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Figure 3. Color-Color selection for redshift z = 4.0 B band 
dropouts. Magenta, red, and cyan X s are model Lyman Break 
Galaxies with different UV continuum slopes at redshifts indi- 
cated along the track. The black bounding box shows the LBG 
selection criteria in color space. Density contours are DLS galaxies 
with a B band S/N > 2 detection. Green points are stars detected 
in DLS. The black track is the evolution of an Eo galaxy CWW 
template from redshift < z < 2. The blue track is a sample 
quasar spectrum redshifted between z = 3 and z = 5. 



distribution of the individual bins, where we have weighted 
each p(z) in a given bin by its corresponding Pi value. Values 
relevant to this analysis for each bin can be found in Table [T] 
To increase the signal to noise of the lensing measure- 
ment, we also optimally weight each foreground galaxy by 
its expected lensing efficiency via 



DiDis 



' . *Da(0, z)Da(z,zlbg) , 

P(Z) - , : -dz 

Py ' D A (0,z LBG ) 



(15) 



where (DiDi s /D s ) is the average, geometric efficiency, p(z) 
is the redshift probability density for a given galaxy, and 
Zmax is the maximum redshift of the p(z) which, given the 
run of BPZ, is z max — 5.0. Da(zi, Z?) is the angular diame- 
ter distance from redshift zi to redshift Z2, and zlbg is the 
redshift of the background LBG sample which is z ~ 4. The 
weight of each foreground galaxy in the correlation is then 
Pi{DiD u /D s ). 



3.3 Lyman Break Galaxy Selection 

To cleanly select a sample of high redshift source galax- 
ies and avoid contamination with the foreground lenses, 
we employ the Lyman B reak criterion as outlin ed in 
iGuhathakurta et all (|l990h and ISteidel et all (| 19991 '). We 
select B band dropout galaxies at 2 ~ 4.0 to use 
as lensed, source galaxies. These gal axies have well un 

(ISteidel et 



derstood luminosity functi o ns (LFs) 



"all 

rthf 



1999 



2007 



Sawicki fc Thompsonl |2006| ; H ouwens fc Illingwort 
van der Burg et al.| [2010). The selection for these galaxies 
in DLS is shown in Figure [3] Using the tracks in color-color 
space for high redshift galaxies with different UV continuum 
slopes, we select a region that avoids low redshift contam- 
inants and stars. The B — V cut is selected to avoid dusty 
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457 K 


0.71 


653 


0.74 



a in Mpc/h 

red cluster galaxies that can mimic the color of LBGs while 
the other cuts avoid both the CWW Elliptical track as well 
as dwarf stars. This color-selection is 

(B - V) > 2 n (V - R) < 2 n (B - V) > {V - R) + 1 (16) 

for z = 4.0 B band dropouts. To ensure that the galaxy has 
dropped out of the B band and is not just due to magnitude 
scatter or deeper/shallower data, we implement a signal to 
noise cut as well as a V band brightness cut that is equal to 
the average depth of B band. 

(S/N) B < 1 n (S/N) v >4n(S/N) R > 5 

nm v <25.5 (17) 

This selection yields ~ 12, 000 LBGs over the whole survey. 

It should be stated that this selection of LBGs is de- 
generate with a selection of quasars (QSOs) at the same 
redshift as evidenced by the QSO track in Figure [3] How- 
ever, as shown by the me asured QSO function within the 
DLS l|Glikman et al-lfeoill ) we expect the number counts of 
QSOs to be two orders of magnitude below that of the LBGs 
for these magnitudes and thus not a significant contaminant. 

For the measurements in >14.2I and >14.3I we impose a 
cut of rriR < 24.8 on the LBGs to ensure they are above 
the LBG completeness limit. However, in we use the 
full range of detected LBG magnitudes to show consistency 
with predicted LBG LFs. 



4 RESULTS 

To esti mate the magnifi c ation -induced cross-correlation we 
use the lLandv fc Szalavl (|l993T) estimator 



(did s ) - (diR s ) - (Rid s ) + (RiR s ) 
(RiR s ) 



(18) 



where {did„) (here we use lower case d to differentiate be- 
tween this quantity and the angular diameter distance, D) is 
the number of pairs between the lensing and source galaxies 
in a given angular bin and Ri is a random sample generated 
using the density and spatial extent of either the lensing or 
source samples. For most of the calculations, the data and 
random points are weighted by some factor Wi making (did s ) 
the sum of the products of the weights for each pair with the 
same set of weights are applied to the random realizations. 
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Figure 4. Lens-source cross correlation due to magnification of background LGBs as a function of LBG magnitude. The left panel shows 
the full angular cross-correlation for each sample against foreground galaxies with rriR < 24 and photometric redshift 0.4 < z p < 1.0. 
In the right panel, we plot the lens-sour ce cross-correlation for a single bin 0.001 < 6 < 0.01 deg for e ach magnitude bin and com pare 
it to the expected scaling for LFs from I Bouwens fc Illingworti] l|2007l l. ISawicki fc Thompson! jjfloj ), Ivan der Burg et al ] (20TCh . and 
ISteidel et ail jl999h . 



Starting with the basic DLS survey footprint, we fur- 
ther mask out regions with bad seeing, high dust extinction 
and shallow depth to ensure a spatially uniform selection of 
foreground and background objects (see Appendix[C]for de- 
tails). This reduces the effective area of the DLS from 20 to 
13.5 deg 2 . We also measure the magnification in each of the 
45 subfields individually rather considering the survey as a 
whole to minimize the effects of varying depth between the 
subfields. This yields a natural, consistent angular scale for 
assembling jack-knife subsamples for the purposes of calcu- 
lating measurement errors. 



4.1 Lyman Break Galaxy Magnitude Bins 

We measure the foreground galaxy - LBG angular angular 
cross-correlation as a function of LBG magnitude. As shown 
in Eq. [6l the magnification signal is expected to scale as the 
slope of the LF: wu ~ a(m) — 1. We split the LBG sample 
into 5 magnitude bins, spanning the range 23 < m^j < 27 
and sampling different parts of the LF with roughly equal 
numbers of LBGs in each bin. The measured angular cor- 
relations Wi s (m) are shown in the left panel of Figure [4] 
The observed density excess of LBGs around foreground 
galaxies is seen to decrease for fainter galaxies, which cor- 
respond to the shallowed part of the luminosity function. 
To explore this more quantitatively, we consider a single 
angular bin running from 0.001 < 8 < 0.01 deg for each 
of the samples and we show the corresponding LBG den- 
sity change in the right panel of Figure [4] We can com- 
pare these values, up to a multiplicative scaling factor, 
with estimates of the slope o f the L B G luminosity function 
from Bouwens fc Illingworth (120071 ) , I Sawicki fc Thompson! 
l|2006h . lvan der Burg et all (|2010h and lSteidel et all ljl999T l. 
As can be seen in the figure, the magnitude dependence of 
the magnification signal is properly recovered. 



— SDSS (M 

• DLS 0.4<z<1.0 




r p [Mpc/h] 



Figure 5. Optimally-weighted magnification-reconstructed mass 
density S as a function of projected radius r p , using all LBG 
candidates from 23 < mjj < 24.8 and foreground lenses between 
0.4 < z p < 1.0. This mass profile is consistent wit h the similarly 
derive d profile from SDSS-based magnification in IMenard et al.l 
l|2010h . 



4.2 Galaxy-Mass Correlation 

In order to combine the measurements from different magni- 
tude bins and maximize the overall S/N of the magnification 
measurement, we apply the o ptimal estimator described in 
IMenard fc Bartelmannl i|2002l ). weight ing each LBG by the 
corres ponding aim) — 1 value from the I Sawicki fc Thompson] 
(2006) LF, which best matched the observed scaling with 
LBG magnitude from Figure [3] Other choices of luminos- 
ity function for the optimal estimator give consistent results 
within the error bars of the measurements, however, they 
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Figure 6. The left panel shows the mass profile reconstruction for each of the 7 photometric redshift bins. For each bin, we calculate 
the tomographic ratio H(z) over the range 0.4 < r p < 1.25 Mpc//i as given in Equation ll3l The right panel shows the covariance matrix 
for these 7 bins. Since TZ(z) is a function of both the lens-source cross-correlation and the lens autocorrelation, any correlation between 
bins in H(z) is a strong indicator of the degree of overlap between redshift bins. We find good photometric redshift bin segregation, as 
required for tomography. 



do tend to bias the measurements low or high depending on 
the relative value of {a — 1) for the LF. From these measure- 
ments we can infer the mean surface mass density around 
the foregrou nd galaxies and com pare the results to those 
obtained by iMenard et al.l (|20T0t ) who used magnification 
measurements of quasars in the SDSS. To convert the mea- 
sured density change into a surface mass density, we first 
compute the average E cr i t for each lens sample based on 
the average lensing kernel for each of the N galaxies in a 
given redshift bin: 



DiDu 



(19) 



With this in hand the conversion from magnification angular 
correlation to surface mass density is 



£(r„) 



wis(r p ) 



(20) 



For the LBG sample of ran < 24.8 and all foreground 
galaxies between 0.4 < z < 1.0, we find a detection S/N 
of ~ 20, as seen in Figure [S] For comparison, we also show 
the results from the mass profile found fo r SDSS galaxies 
via magnification from lMenard et all (|201p| ). The DLS mea- 
surements are at higher redshift (z ~ 0.4 for the SDSS sam- 
ple) and redshift-selected rather than magnitude-selected, 
so some level of disagreement is expected. However, from 
a qualitative standpoint, there is good consistency between 
the two measurements over a large range in scale. 



4.3 Tomography 

We now explore the redshift-dependence of the lensing sig- 
nal. We estimate the galaxy-mass correlation function for 7 
subsamples of foreground galaxies with varying redshifts as 
described in £13.21 The left panel of Figure[6]shows the values 
of £(r) as a function of lens redshift. The observed redshift 



trend is due to a number of effects: due to our brightness 
and color selection, galaxies with different types and masses 
are selected in different redshift bins. In addition the shape 
of the redshift distribution differs from a bin to another one. 

As mentioned previously, we can combine the galaxy- 
mass cross-correlation and galaxy auto-correlations through 
the quantity 7Z(z) introduced in N2.ll to obtain a quantity 
which, on large scales, does not depend on the galaxy prop- 
erties. To measure 7Z(z), we use a single bin of constant 
physical radius over 0.4 < r p < 1.25 Mpc/h for each pho- 
tometric redshift bin. This choice of spatial scales is limited 
by the range over which 7Z(z) is independent of galaxy bias 
(see Appendix|X]|, as well as the acceptable scales where Wi s 
and wu have good S/N. 

Before discussing the TZ(z) measurements themselves, 
we consider the normalized correlation matrix for TZ(z) in 
the right panel of Figure [|J] Since the covariance between 
redshift bins in TZ(z) is a combination of the covariance for 
both the magnification and foreground autocorrelation sig- 
nals, the level of correlation between redshift bin pairs is ex- 
pected to be sensitive to the redshift overlap between those 
bins. The correlation which we find is consistent with what 
we expect from the distributions in Figure [2] This suggests 
that we have the level of photometric redshift segregation re- 
quired to observe tomographic lensing. We also find that this 
is robust under a variety of cuts in P-threshold (see i )3.2[) 
and photometric and radial binning. In using this matrix 
for deriving co nfidences, we c orrect the inverse covariance 
as described in lHartlap et aD l|2007n . 

We present measurements of 1Z{z) in Figure [7] for each 
of the 7 photometric redshift bins, along with a number of 
theoretical models. The left panel compares the observed 
7t(z) with several flat cosmologies with increasing S2 m and 
fixed as and h. For a fixed set of redshift bins, the amplitude 
of IZ(z) scales quite strongly with Q m since the overall am- 
plitude of the lensing signal is directly proportional to the 
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Figure 7. Tomographic ratio (TZ(z)) for photometric redshift bins over 0.4 < z < 1.0 on scales 0.4 < r p < 1.25 Mpc/h. The left panel 
shows H(z) for cosmologies over a range of f2 m values. The shaded region shows a variation of ±0.1 with respect to the fiducial value of 
a g . For fixed ag, we reject a matter dominated universe at > 7. 5a and prefer a universe with f2 m = 0.3. The right panel compares the 
measured TZ(z) with two potential null models. We reject the scrambled photometric redshift model at 3.8a, however the measurements 
disfavor the no-tomography theoretical null model at only 2.6a. 



mean matter density. We also plot a shaded region around 
fim = 0.3, varying ag by ±0.1 around the feducial value of 

0. 8 to give the reader a sense of the usual erg, Q, m degen- 
eracy. We find that the measurements reject a flat, matter 
dominated model at > 7.5a, preferring a cosmology with 

~ 0.3 (x 2 /v = 1.7) for a fixed value of og. 

However, this consistency does not necessarily verify 
that we are observing tomographic lensing. To do so, we 
consider two potential null models, one observational and 
one theoretical. The observational null model assumes that 
the photometric redshift binning has failed completely and 
the redshift bins are essentially random subsamples of the 
full foreground lens sample. The theoretical model assumes 
that K(x) from Equation[S]has been replaced by a constant, 

1. e. that the lensing kernel is flat as a function of redshift. 
There is no obvious value for us to use; therefore we set this 
constant to the mean kernel value of the lensing sample. 

We plot the results of these tomographic null tests in 
the right panel of Figure [7] The measurements reject the 
scrambled lens redshift null test at 3.8ct and prefer a ACDM 
universe at 4.8a, as one would expect based on the correla- 
tion matrix in Figure O We reject the fixed constant theo- 
retical model at 2.6a and prefer the concordance cosmology 
to this model at 3. 0<t. This mild rejection of the theoretical 
null model is not entirely surprising given that the lensing 
kernel across these samples is relatively flat, as seen in Ta- 
ble [1] Had we been able to make measurements at either 
higher or lower redshift, rejecting the null theoretical model 
might have been possible. We explored pushing beyond the 
lens redshift limits shown in the right panel of Figure [2l but 
the combination of small survey angular extent and photo- 
metric redshift degeneracies made measurements at those 
redshifts unreliable. 



5 DISCUSSION AND CONCLUSIONS 

Using data from the Deep Lens Survey, we have shown a 
robust, S/N > 20 detection of the cosmic magnification 
of background LBGs by foreground galaxies , producing a 
halo m ass profile similar to that observed in lMenard et al.l 
(2010). In addition, we used photometric redshifts to divide 
the foreground sample into 7 redshift-selected bins, recover- 
ing a strong detection in each as well as demonstrating that 
the measurements had the expected level of redshift segre- 
gation. By combining these lensing measurements with the 
measured angular autocorrelation in each foreground sam- 
ple, we effectively de- biased the measurement, isolating the 
cosmological signal 7t(z). 

We tested these measurements of 1Z in two ways. First, 
we compared it to the expected signal for flat universes with 
fixed h and a g and varying Q m . In this test, we found that 
the results were consistent with the concordance cosmology 
and rejected a flat matter dominated universe at > 7.5ct. 
Second, we tested the tomographic nature of the TZ(z) mea- 
surement by considering two null models: an observational 
null model where we assumed that, despite indications oth- 
erwise, the photometric redshifts failed to separate the fore- 
ground samples; and a theoretical null model where we as- 
sumed that the lensing kernel was a constant set to the mean 
value of the kernel for this sample. We reject the observa- 
tional null model at 3.8ct but the theoretical null model was 
only rejected by the data at 2.6ct, owing to the fact that the 
lensing kernel only varied by ~ 10% over the 7 redshift bins 
and higher and lower redshift samples were unreliable due 
to the limitations of the survey. 

From these results we can draw several conclusions: 
First, that magnification can yield high signal to noise with- 
out the necessity for complex shape measurements; second, 
that magnification can be combined with autocorrelation to 
suppress galaxy bias and constrain the integrated dark mat- 
ter distribution. While the method in this paper obviates 
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the need for it, fitting a full halo model to the data is fea- 
sible and would allow this measurement to both constrain 
cosmology directly as well as measure the galaxy bias. Up- 
coming papers will address this extension. 

Future work will include combining with measurements 
of galaxy-galaxy lensing measured using the same fore- 
ground lens sample. This will leverage the additional con- 
straining power of both the galaxy bias and the multi-sourc e 
and multi-lens weak lens tomography (jvan W acrbckc 2010h . 
With the large range of redshifts as well as the ability of the 
magnification to utilize unresolved objects, we can probe 
mass-luminosity relations at z ~ 1.0 as well as the d ust con- 
tained in high redshift galaxies (|Menard et al.ll2010T l. 
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the halo, M, is above some threshold mass Mo. For the satel- 
lite galaxies, the expected number of galaxies in a halo of 
mass M is 



(N a \M) oc 



M if M < 3M 



M if M > 3M 
The average number of total galaxies in a halo is then 
(N\M) = {N C \M} + (N S \M) 



(Al) 



(A2) 



For values of Mo we use the measured values in three differ- 
ent luminosity samples at z ~ 0.5 as observed in the DLS 
by galaxy-galaxy lensing (private communication A. Choi). 

The Z07 HOD also uses a central/satellite galaxy model 
as a function of halo mass. The functional form of the central 
term is 



(N C \M) 



1 + erf 



log(M) - log(M min 

ClogM 



(A3) 



where M m i n is the minimum mass for a halo to have one 
galaxy and o"i og M is the width of the central galaxy turn on. 
The satellite galaxy term is 



(N S \M) = 



M -Mp 
M[ 



(A4) 



where Mo is the minimum mass for a halo to host satellite 
galaxies (note this is distinct from the Mo in the Mandel- 
baum model), and M[ is the mass differential at which a 
halo is expected to have one satellite galaxy. The average 
number of galaxies occupying a halo of a given mass is then 



(N\M) = (Nc\M)(l + (N s \M)) 



(A5) 



We use the measu red parameters for this HOD model from 
IWake et al.l (|201 lh for galaxies at z — 1.1 in three bins of 
stellar mass. 

The form of t he HOD second moment we use is from 
IZheng et al.l (|2005h . 



(N(N - l)\M) = 2(N B \M) + (N S \M) 



(A6) 



APPENDIX A: TOMOGRAPHIC RATIO BIAS 
DEPENDENCE 

To estimate the bias dependence o f 1Z we imp lement a 
halo model pow e r spe ctrum code l|Seliakl |2000| ) with a 
ISheth fc Tormenl (2002) mass function. Within this model, 
we use two different halo o ccupation distributions (H ODs), 
based on forms given by iMandelbaum et al l (|2005l ) and 
I Zheng et all (|2007m (M05 and Z07, hereafter, respectively). 
The galaxy bias for the magnification and autocorrelation 
are determined by the first {(N)) and second moments 
({N(N - 1))) of the HOD respectively. So long as the ra- 
tio between these two moments is Poissonian (< TV > 2 ~< 
N(N — 1) >), then any scale dependence in the bias for the 
samples will cancel out, leaving Equation [13] bias indepen- 
dent. 

The M05 HOD breaks the expected number of galax- 
ies per halo into satellite ((7V S |M)) and central ({N C \M)) 
galaxies. A halo has a single central galaxy if the mass of 



where N„ is the expected number of satellite galaxies in a 
halo of mass M. 

Table IA1I presents the parameters for the HODs used, 
along with the mean redshift and linear galaxy bias. As dis- 
cussed in §3] the redshifts for the lensing samples are brack- 
eted by the galaxie s used in t he DLS galaxy-gala xy lensing 
analysis Ichoi et all (|2012T ) and lWake et alj l|201lf ), so we ex- 
pect the results to be bounded by these models. Since the 
galaxies are selected by photometric redshift within fixed 
apparent magnitude bounds, we expect them to be slightly 
biased relative to field galaxies, with the bias increasing as 
we go to higher redshifts. As shown in Figure [l] the various 
HOD models track the overall shape of the dark matter- 
only model reasonably well and are within ~ 20% of the 
dark matter ratio for all of the redshift bins even at scales 
as small as 0.5 Mpc/ft. This suggests that Equation 1131 is 
adequately galaxy bias independent on these small scales 
and, to within the precision of this measurement, probes 
the magnification due to dark matter directly. 
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Table Al. Mean Redshift and HOD Bias 



HOD Model Name 


(z) 


M 


b 9 


DLS A 


0.58 


0.13 


1.60 


DLS G-G Lensing DLS B 


0.58 


0.02 


1.43 


DLS C 


0.53 


0.01 


1.27 



HOD Model Name 






M[ 


bg 


Wake A 


1.1 


0.31 


1.45 


2.59 


Wake et al. (20U) a Wake B 


1.1 


0.17 


0.76 


2.37 


Wake C 


1.1 


0.10 


0.42 


2.17 



All Masses are in 10 13 Mq /h. 
"Several HOD parameters are fixed for this measurement, they are: Mq = M m i n ; <xi og M = 0.15; a = 1.0. 



APPENDIX B: DEEP LENS SURVEY 
CALIBRATION & PHOTOMETRIC REDSHIFTS 

Bl Photometric Calibration 

The DLS image processing and photometry are described 
in detail in an upcoming data release paper (Wittman 
et al. in prep.); here, we give a brief overview. The im- 
ages were de-biased, flat-fielded, and sky-subtra ct ed us - 
ing the procedures described in IWittman etafl (|2006h . 
Next, we used t he global, linear least-squares algorithm 
know n as ubercal l|Padmanabhan et al ] |2008l ; IWittman" et al.l 
l201ll ) to correct for residual fiat-fielding errors, after which 
the corrected images were cross-registered and photometri- 
cally stacked. Object detection was done with the R band 
(the deepest and highes t-resolution filter) using SExtractor 
l)Bertin fc Arnoutslll996t ). with magnitudes measured in the 
other bands regardless of S/N in that ban d. We used a modi- 
fied version of ColorPro (|Coe et alj|200rjfr to provide robust- 
ness against variations in seeing because photometric red- 
shifts require matched aperture photometry. To measure a 
given color (e.g. B — R), ColorPro convolves the R band 
image to match the B image and then uses matched aper- 
tures on the seeing-matched images. The resulting photom- 
etry showed spatial variations of up to ~ 0.04 mag in 2 band 
and ~ 0.03 mag in B and V bands, based on the stellar lo- 
cus in color-color space. We corrected for this by applying 
zero-point shifts to each subfield to make the stellar locus 
consistent across subfields and across fields. This was gener- 
ally a straightforward correction, but the varying shape of 
the stellar loci with galactic latitude and longitude compli- 
cated the overall shifts of fields Fl and F2 with respect to 
each other and the rest of the survey. To determine these 
shifts, we used overlapping p arts of the Sloan Digital Sky 
Survey (|Eisenstein et al]|201ll ) as the basis for common cal- 
ibration and removed the resulting shifts, which ranged from 
0.00 to 0.05 mag depending on field and filter. We estimate 
that the remaining spatial variation in the photometry is at 
the level of 0.02 mag or less. 

B2 Photometric Redshifts 

The photometric redshifts are based on the BVRz photom - 
etry obtained from the ColorPro and BPZ i|Bemtezl I2OO0T ) 
software packages. We replaced the standard templates with 
a set optimized in a method similar to that described in 



lllbert et al.l (|2006t ). By using spectroscopic samples from the 
SHELS survey ( private communic ation M. Geller) and the 
PRIMUS survey (|Coil et all (|201lh . Cool et al. in prep.) that 
overlap the DLS footprint, we divided the galaxies into six 
galaxy types [Elliptical, Sbc, Scd, etc.) and then adjusted 
the SED templates to match the median rest-frame fluxes 
observed in the DLS photometric data as a function of wave- 
length. This procedure matched the colors of galaxies to the 
observed data, reducing "template mismatch" bias and vari- 
ance. 

We also employed a modified version of the type-redshift 
prior used in BPZ. Beginning with the spectroscopic data 
from SHELS, we fit the prior to the observe d P(z\T, m) dis- 
tribution using the form from lBenite3 (|2000T ). This prior was 
extended to f ainter magnitud es using the VVDS spectro- 
scopic sample (|Le FevrAooBT ). A more detailed description 
of the DLS photometry and photometric redshifts will be 
given in Thorman & Schmidt (in prep.). 

APPENDIX C: SOURCES OF SYSTEMATIC 
ERRORS 

CI Observational Contaminants 

Since both the foreground and background samples are 
drawn from the same underlying imagery, inhomogeneities 
in the detection and selection process due to variations in 
the seeing, RMS sky noise, and extinction can induce cor- 
relations between the two samples which are unrelated to 
the lensing magnification signal we are trying to detect. To 
reduce the effect of these biases, we exclude regions around 
bright stars and CCD bleeds. We also mask out the regions 
at the edges of each subfield where the number of stacked 
exposures can be as low as ~ 4, compared to ~ 20 at the 
center of a subfield. This brings the nominal area down from 
a full survey area of 20 deg 2 to 13.5 deg 2 . 

To estimate the contribution of the various sources 
of systematic error to the angular correlations, we cross- 
correlate both the foreground lenses and background LBG 
against the seeing, RMS sky noise, and extinction for each of 
the 4 bands. These correlations are plotted in Figure [ClI for 
the R band. Other bands have similar values for the correla- 
tion and are not shown here. Because we undertake a cross- 
correlation between a background and foreground sample, 
we can tolerate a certain degree of systematic correlation in 
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Figure CI. Angular correlations between the foreground lens (left) and LBG (right) samples against various potential sources of 
systematic errors, th e RMS sky noise (Sky), the mean seeing (See) full width half maximum, and extinction (Ext) from galactic dust 
jSchlegel et al.lll99Sh . We compare these systematics with the respective sample autocorrelation functions and the magnification signal. 



one sample, so long as the other is uncontaminated. From 
the cross-correlations plotted in Figure ICll we see that this 
is the case for the foreground lenses against the systematic 
sources as they do not show spatial structure, are gener- 
ally consistent with zero, and are small compared to the 
amplitude of the magnification and autocorrelations. Since 
the foreground systematics are uncorrelated and the LBG 
systematic cross-correlations are small, the effect of these 
systematic errors on these measurements is small. 



C2 Redshift Contaminants 

As stated in iJSl the most dangerous contaminant in the 
LBG sample are low redshift, dusty, red galaxies. If the 
LBG sample is contaminated by these low redshift galax- 
ies at even 10%, the signal could be due to physical clus- 
tering and not magnification. To test for contamination at 
low redshi ft we cross-correla te the LBG sample against the 
PRIMUS (|Coil et all <|201lf) . Cool et al. in prep.) spectro- 
scopic sample that overlaps the DLS footprint, in a single, 
large radial bin (100 kpc to 1 Mpc). The PRIMUS sample 
has ~ 8000 spectra with roughly ~ 2500 LBGs in the in- 
tersection of the PRIMUS footprint with the 13.5 deg 2 used 
in thi s anal ysis. This method is similar t o that of iNewmanl 
l|200Sl ) and' lMatthews fc Newman! (|2012l ) and exploits the 
physical clustering of galaxies to recover redshift distribu- 
tions. While this spectroscopic sample does not allow us to 
recover the redshift distributions beyond z ~ 1, we can uti- 
lize it as a null test for contamination from low redshift 
galaxies, (see Menard et al. 2012 in prep, for full details of 
the algorithm.). The PRIMUS galaxies induce a correlation 
on the LBGs through magnification, however, this correla- 
tion will nulled out when plotted as an over-density relative 
to the average density over the whole redshift range as it 
varies slowly (See Table [TJ and will therefore not have the 
same signature as the clustering correlation in this test. 

If we are truly selecting Lyman Break Galaxies at red- 
shift 2 ~ 4 then we should see no significant cross-correlation 




-o.a - 



Figure C2. Test for low redshift galaxy contamination of 
the LBG sample , using the PRIMUS spectroscopic rcdshifts 
dCoil et al.l <201lD . Cool et al. in prep.). Plotted is the fractional 
over-density of LBGs in a 100-1000 kpc radial bin around the 
z < 1 spectroscopic objects. If there were a significant number of 
low redshift galaxy contaminants of the LBG sample, this corre- 
lation would show some structure with redshift. As it does not, 
we can be confident that the LBG sample is free of low redshift 
contaminants. The induced correlation by magnification of the 
LBGs by the PRIMUS galaxies will be mostly nulled out as the 
over-density is computed with respect to the average density over 
the redshift range plotted and the amount of magnification over 
this range varies slowly (See Table [TJ. 



between these low redshift spectroscopic samples and the 
candidate LBGs. Figure [C2l shows the mean over-density of 
LBGs in an annulus surrounding each spectroscopic object 
in PRIMUS as a function of redshift (error bars are from 
a spatial jackknife). The lack of coherent structure above 
that expected from an uncorrelated background sample im- 
plies that there is no detectable contamination of the LBG 
sample by low redshift galaxies. 
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